/3-Mn: Emergent Simplicity in a Complex Structure 
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We investigate low-temperature spin correlations in the metallic frustrated magnet /3-MnCo. 
Single-crystal polarised-neutron scattering experiments reveal the persistence of highly-structured 
magnetic diffuse scattering and the absence of periodic magnetic order to T = 0.05 K. We employ 
reverse Monte Carlo refinements and mean-field theory simulations to construct a simple effective 
Hamiltonian which accounts for the magnetic scattering. The interactions we identify describe an 
emergent spin structure which mimics the triangular lattice antiferromagnet. The observation of a 
simple collective magnetic state in a complicated crystal structure is surprising because it reverses 
the established paradigm of elaborate emergent states arising from many-body interactions on sim- 
ple lattices. We suggest that structural complexity may provide a route to realising new states of 
correlated quantum matter. 
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In frustrated magnets, geometry and magnetic interac- 
tions can conspire to suppress the formation of long-range 
magnetic order. A theme common to the resulting dis- 
ordered states is emergent behaviour: topical examples 
include quasiparticles resembling magnetic monopoles in 
spin ices [1], composite spin-orbital degrees of freedom 
[2], and spinon excitations in quantum spin liquids [3]. 
Typically, investigations of collective behaviour have fo- 
cused on insulating materials in which magnetic ions oc- 
cupy a frustrated lattice such as kagome or triangular. 
In such cases, remarkably complex emergent states can 
arise from a simple geometry and set of magnetic inter- 
actions. In this Letter, we demonstrate a reversal of this 
paradigm: the presence of simple emergent spin struc- 
tures in Co-doped /3-Mn, a metal which, although chem- 
ically trivial, has a crystal structure well-known for its 
complexity 0]. 

Pure /3-Mn is a unique example of an elemental quan- 
tum spin liquid, which exhibits a host of largely unex- 
plained magnetic phenomena, including geometrical frus- 
tration large antiferromagnetic spin fluctuations [6], 
and non- Fermi liquid behaviour [7^. However, all previ- 
ous experimental studies of pure /3-Mn have been limited 
by the lack of single-crystal samples. Here, we investi- 
gate two large single crystals of /3-Mn, doped with 12% 
and 20% Co. The availability of single-crystal samples 
means that /3-MnCo is a much better candidate for ex- 
periments than elemental /3-Mn. Yet element and alloy 
are in many respects similar: the lattice parameter re- 
mains almost unchanged on Co-doping, while spin fluc- 
tuations are damped in the alloy but do not freeze out 
[SI [9]. The /3-Mn structure is primitive cubic (space group 
P4i32), with a unit cell containing 20 Mn atoms divided 



over two inequivalent crystallographic sites (8c and 12d) 
[3]. Neutron diffraction reveals that Co substitutes al- 
most entirely onto the 8c site (TU]. We will show that 
all our data can be explained in a model in which the 
8c site is non-magnetic, as has already been proposed 
for elemental /3-Mn [SI [TT]. Accordingly, in what fol- 
lows we consider the magnetism of /3-MnCo in terms of 
the \2d site only. Connections between neighbouring Mn 
atoms form a three-dimensional "hyperkagome" network 
of corner-sharing triangles [Fig. la] — a topology which 
has led several authors to identify the possibility of geo- 
metrical frustration [3 IHl [H] . 

Our neutron scattering experiments provide strong ev- 
idence that frustration indeed plays a central role in /3- 
MnCo. Measurements were performed on the D7 diffrac- 
tometer at the Institut Laue-Langevin, Grenoble, using 
xyz polarisation analysis to isolate the magnetic compo- 
nent to the scattering [Ij]. The sample with 20% Co 
was a 17.2 g single crystal prepared using the Bridg- 
man technique. Measurements of this sample at tem- 
perature T — 1.5 K are shown in four reciprocal-space 
planes in Fig. 2 (top left panels). These data reveal 
highly-structured diffuse scattering patterns — a hallmark 
of frustrated magnetism [14] . Furthermore, tempera- 
ture scans reveal that this diffuse scattering remains un- 
changed over three orders of temperature magnitude be- 
tween T = 0.05 and 50 K, with no sign of the development 
of long-range magnetic order at T = 0.05 K (see SI). Im- 
portantly, measurements on the 12% sample show that 
the form of the diffuse scattering is robust to Co concen- 
tration, indicating that the magnetic short-range order is 
unlikely to result from Co-doping, but is instead an in- 
herent property of /3-Mn (see SI). These results suggest 
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FIG. 1: (a) The /3-Mn crystal structure (space group P4i32), 
showing magnetic Mn atoms (12d site) only. Frustrated 
nearest-neighbour connections are shown as grey or burgundy 
lines. While all connections shown are of identical length to 
within 1% (see SI), only connections shown the same colour 
can be brought into coincidence by symmetry, (b) Local 
Mn environment. The six nearest neighbours of each Mn 
are arranged in a "windmill" geometry, with each triangu- 
lar "blade" perpendicular to a (111) axis. The local x, y, and 
z (out of the page) axes are defined for each Mn atom relative 
to the burgundy-coloured triangle. This choice of coordinates 
aligns the x axis with the two-fold rotation axis which is the 
only point symmetry element of the 12d site. 



two key questions. First, what is the microscopic origin 
of the persistent magnetic disorder? Second, what is the 
nature of the real-space spin correlations which give rise 
to structured diffuse scattering? 

The simplest model which might address these points 
considers only antiferromagnetic interactions between 
nearest-neighbour spins [Fig. lb] [5|. At the mean-field 
level, these interactions lead to a macroscopic degener- 
acy of magnetic ground states f5]. However, the diffuse 
scattering calculated from this model disagrees markedly 
with our data [Fig. 2, bottom left panels]. In order to 
develop a more sophisticated model, we employ a reverse 
Monte Carlo (RMC) approach [15], in which the orien- 
tations of classical spin vectors in a large configuration 
are refined in order to fit experimental data. The sum of 
squared residuals minimised during the RMC refinements 
is given by 



(1) 



where /(Q) is the magnetic scattering intensity at 
reciprocal-space position Q, superscript calc and expt 
denote calculated and experimental values, s is a refined 
intensity scale factor, cr(Q) an experimental uncertainty, 
and W an empirical weighting factor. Subscript d de- 
notes different datasets, namely the four complete single- 
crystal reciprocal-space planes shown in Fig. 2 and the 
powder neutron scattering data of Ref . [S] . The magnetic 
scattering intensity is calculated using the standard ex- 
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FIG. 2: Magnetic neutron scattering data were collected in 
four reciprocal-space planes, clockwise from top left (001)*, 
(110)*, (112)*, and (111)*. In each plane, the top left panel 
shows T = 1.5 K experimental data, the top right panel shows 
the RMC fit to data, the bottom right panel shows the scatter- 
ing for the J1-J5 model described in the text, and the bottom 
left panel shows the scattering for antiferromagnetic Ji inter- 
actions only. Intensity scales are chosen such that the maxi- 
mum value of /(Q) in each plane is equal to 1, and m3m Lane 
symmetry has been applied to experimental data and RMC 
fits. The Ji and J1-J5 models are calculated using a mean- 
field theory at simulation temperature T = LOIT^*^^, where 
is the mean-field transition temperature. Whereas the 
Ji-only model is a very poor description of the data, both 
the RMC fit and the J1-J5 model give good agreement with 
experiment. 
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where /(|Q|) is the Mn^+ magnetic form factor |16| . 

= Si — [(Sj • Q)Q]/|Q|^ is the component of spin 
Si perpendicular to Q, and the position of spin Sj. 
Refinements were performed with spin configurations of 
size 10"^ unit cells (12000 spins) using periodic boundary 
conditions; full details are given as SI. The fit to data ob- 
tained by our RMC refinements is shown in Fig. 2 (top 
right panels). In the refinements we constrain the 8c 
site to be non-magnetic — an assumption supported by 
the excellent agreement obtained between data and fit, 
but which is in apparent conflict with the conclusions 
of a previous RMC study [9]. We have since identified 
computer-programming errors in the RMC code used in 
that previous study, which are entirely responsible for 
the misdiagnosis of a magnetic contribution from the 
8c site (see SI). In the present study, the use of single- 
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crystal samples together with an improved RMC algo- 
rithm enables us to develop a fundamentally more ro- 
bust and detailed description of spin correlations in /3- 
MnCo. So, access to these new RMC spin configura- 
tions finally allows us to calculate the true real-space spin 
correlation functions responsible for the observed diffuse 
scattering. We consider first the radial spin correlation 
function (S(0) • S(r)) shown in Fig. 3a. As anticipated 
from the absence of neutron scattering intensity near 
Q = 0, strongly antiferromagnetic nearest-neighbour 
correlations are present. However, we find that ferro- 
magnetic further- neighbour correlations are prominent, 
and note that ferromagnetic correlations between fifth- 
nearest-neighbours are as strong as nearest-neighbour an- 
tiferromagnetic correlations. This finding is perhaps un- 
expected given conventional descriptions of /3-Mn as a 
frustrated antiferromagnet [Gj, but may nonetheless be 
consistent with spin-polarised density-functional theory 
studies, which predict ferromagnetic correlations within 
the 12 d Mn sublattice [Hj . 

To complement the model-independent RMC 
approach — which identifies the strongest correla- 
tions — we seek now to develop a Hamiltonian for 
/3-MnCo that captures the key magnetic interactions. 
Analysis of the RMC spin configurations indicates that 
the degree of spin anisotropy is small, suggesting that 
an effective Heisenberg model is appropriate (see SI). 
The Heisenberg Hamiltonian can be written 
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H — ^ ^ Jab (R-, R') Sr_q • Sr'^ 



(3) 
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where Jab (R-, R-') = Jn is a coupling between nth. neigh- 
bours, and each classical spin vector Sr ^ = is labelled 
by a site index a € {1, 12} and a primitive lattice transla- 
tion vector R. Here, the Heisenberg Hamiltonian is used 
as a phenomenological model, in which magnetic inter- 
actions between both localised and conduction electrons 
are renormalised into a small number of effective param- 
eters. This effective local-moment approach has been 
widely used to study metalHc systems (see, e.g., JT]). 
We use a mean-field theory [TB] to calculate the param- 
agnetic neutron scattering pattern for a given set of J„. 
The Fourier transforms of the Jab (R, R') are written as 
a 12 X 12 Hermitian matrix J (q), whose eigenvalues A (q) 
and eigenvectors U (q) are related to the magnetic neu- 
tron scattering intensity by the expression 
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Aa(q)/T) 
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where a, a, 6 G {1, 12}, and Q = q -I- G, with q a vector 
in the first Brillouin zone and G a primitive reciprocal 
lattice vector [TS]. The interaction parameters J„ were 
fitted to the data shown in Fig. 2 using a least-squares 
approach (see SI). In order to determine the minimum 
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FIG. 3; (a) Radial spin correlation function (S(0) ■ S(r-)) 
calculated from RMC spin configurations, (b) Variation of 
the sum of squared residuals as a function of the number 
of interaction parameters J„ allowed to vary when fitting a 
Heisenberg model (Eq. 1) to the data. The minimum for 
Ji + Jn (red line) at n = 5 indicates a good fit to data, which 
is not significantly improved by allowing a third parameter to 
vary (green hue) or by including all Jn up to n = 12 (blue 
line). 



number of J„ sufficient to describe the data, we plot the 
dependence of the sum of squared residuals on the 
number of J„ allowed to vary. The results of this proce- 
dure are surprising yet unambiguous [Fig. 3b]: only two 
interaction parameters are necessary to fit the data, with 
the inclusion of up to 10 further parameters leading only 
to marginal additional improvements in the fit. The two 
essential parameters are antiferromagnetic Ji and ferro- 
magnetic J5, in the ratio J5/J1 ~ —0.65. We note that 
these interactions correspond to the largest magnitudes 
of (S(0) • S(r)) in Fig. 3a. The high quality of the fit to 
neutron scattering data obtained [Fig. 2 (bottom right 
^ panels)] is remarkable given that the model involves just 
two parameters. 

This simple J1-J5 Hamiltonian helps explain the ab- 
sence of magnetic order in /3-MnCo. In a conventional 
magnet, the largest eigenvalue of the Fourier transform 
of the interaction matrix with elements Jab(R, R') — 
denoted A„iax(q) — bas a global maximum at some 
wavevector q = qord in the first Brillouin zone. At the 
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mean-field level, this is the propagation vector of the first 
ordered state [19]. By contrast, in a frustrated magnet, 
Amax(q) is flat in q, preventing the formation of an or- 
dered state [20] • In real materials, a completely flat spec- 
trum of Ainax(q) is never obtained. However, the spec- 
trum of Aniax(q) in /3-MnCo is nearly flat across a large 
surface in q-space (full details are given as SI) , suggesting 
that frustration of the J1-J5 Hamiltonian is a key factor 
contributing to the lack of order. 

Having obtained a credible interaction model, we pro- 
ceed to determine how these interactions are manifest in 
the spin correlations obtained from RMC fitting. The ra- 
dial spin correlation function (S(0) • S(r)) shown in Fig. 3 
does not suffice for this purpose, because the radial av- 
erage obscures the relationship of the crystal structure 
to the spin correlations. We therefore consider the three- 
dimensional spin correlation function (S(0) • S(r)) . When 
moving from one spin to the next in this configurational 
average, we take account of the crystal symmetry by 
defining local coordinates (x, y, z) which reflect the crys- 
tallographic equivalence of each Mn atom [Fig. 1]. The 
relationship between Ji and J5 interactions is shown in 
Fig. 4. An arbitrary central Mn atom is coupled by fer- 
romagnetic J5 interactions along helical chains [Fig. 4a]. 
Representing these chains as rods oriented along the heli- 
cal axis maps a set of helices onto a triangular rod lattice. 
In this picture, the entire structure is described as four 
such rod lattices, oriented along the cubic (111) direc- 
tions [3] . The central Mn is coupled by antiferromagnetic 
Ji to two atoms within the same rod sublattice; its four 
other nearest neighbours belong to the other three rod 
sublattices [Fig. 4c]. 

The mapping of the /3-Mn structure onto a rod packing 
suggests considering the spin correlations in two different 
orientations: along the J5 rods, and in the plane perpen- 
dicular to them. We consider first the spin correlations 
along J5 rods, denoted (S(0) • S(r^)). This represents 
the correlation of our reference Mn spin with its neigh- 
bours at distance within the central rod in Fig. 4a. 
Values of (S(0) • S(rz)) are shown in Fig. 4b: they are 
always ferromagnetic, with a simple exponential decay as 
a function of r^. This trivial dependence means that 
it is possible to interpret a rod as behaving like a sin- 
gle collective object, and to calculate the summed spin 
correlation value (^(0) ' ^(''z)) ^'^^ ea,ch rod within 
the rod packing oriented along z. These sums do not di- 
verge because the ferromagnetic correlations have finite 
length. Values of J2r i^i^) ' ^(''2)) are shown in Fig. 4d 
(right panel). The sign of X]r-(^('^) ' ^('^z)) alternates 
throughout the triangular rod packing, because adjacent 
rods are coupled by Ji. This pattern strongly resembles 
the (S(0) • S(r)) calculated for a simple antiferromagnetic 
Heisenberg model on the simple triangular lattice [21]l22j. 
shown in Fig. 4d (left panel). We therefore suggest that 
the comparable strength of Ji and J5 leads, over short 
length scales, to frustration of multi-spin rods on a tri- 
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FIG. 4: (a) Ferromagnetic J5 interactions couple spins along 
helical chains, shown as thin red lines. The direction of propa- 
gation of each chain is shown as a thick red rod. The complete 
crystal structure can be described as four rod sublattices, 
each of which describes a triangular rod lattice. For clarity, 
only a single rod sublattice is shown, (b) Spin correlations 
(S(0) • S{rz)) along J5 rods are always ferromagnetic, and 
decay exponentially with distance along the rod with cor- 
relation length ~5 A. (c) Antiferromagnetic Ji interactions, 
shown as thin lines, couple spins from adjacent rods. An ar- 
bitrary central Mn atom, shown as a large red sphere, is cou- 
pled to two atoms within the same rod sublattice (coloured 
red); its four other nearest neighbours belong to the other 
three rod sublattices (coloured green, blue and yellow), (d) 
The right panel shows spin correlations {S{0)-S{rxy)) perpen- 
dicular to J5 rods. Values shown represent a sum along each 
rod, as described in the text. The spin correlations alternate 
between antiferromagnetic and ferromagnetic throughout the 
triangular rod lattice. The left panel shows spin correlations 
(S(0) ■ S(r)) for a frustrated antiferromagnetic Heisenberg 
model on the simple triangular lattice (see SI). 



angular lattice. This emergent behaviour contrasts with 
the frustration of single spins on a hyperkagome lattice 
which would result from a Ji-only model [S]. 

The identification of strong magnetic frustration and 
emergent spin structures in such a chemically-simple sys- 
tem provides a valuable experimental reference point 
against which to benchmark developments in the theory 
of unconventional metals [331 IM] ■ Our results may pro- 
vide a tantalising glimpse at the behaviour of pure /3-Mn: 
the magnetic diffuse scattering is robust to Co concen- 
tration in the range for which single crystals can be pre- 
pared, while the diffuse scattering from powder samples 
remains essentially unchanged for the entire range of Co 
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concentration from zero to 20%, with only an increase in 
the static magnetic moment evident on doping [6l[9]. But 
perhaps our key result is to show how typical behaviour — 
the emergence of complex behaviour from simple struc- 
tures, as has previously been observed in metallic systems 
[551 [2S] — can be reversed in a structurally-complex sys- 
tem. There are hints here of a more general phenomenon: 
the emergence of simple collective states within complex 
networks has also been observed in fields as diverse as 
neural signalling [37] and the dynamics of amorphous 
materials [23]. Hence a focus on structural complexity 
may prove an important experimental strategy for real- 
ising many of the as-yet unobserved states of correlated 
quantum matter |29j . 
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